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I. Introduction 


Written in generalized curvilinear coordinates £ — £(x, y), P — vi-^i y) fhe inviscid 
Navier-Stokes equations are 

■Q t Q + •O^F + = 0 (1) 
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where, p, u, v,p, e are conventional physical properties, and 

U = (PI + tyV, V = T] X U + T] y V. 

Compared to the Navier-Stokes equations written in Cartesian coordinates (x,y), except 
with the additional inclusion of transformation Jacobian J, (U,V) are the replacements 
of velocity (u, v), and hence are called “contravariant velocity”, or “contravariant velocity 
components” based on their associated with the contravariant vector. 

The term “contravariant velocity component” is used uniformly in computational 
fluid dynamics (CFD) community and publications, in general without further statements 
or explanations, for instances, Refs. 1-2. In Ref. 3, while applied for body surface, it 
stated that 


“ the contravariant velocities are the decomposition of velocity vector ^ into 

components along the £ coordinate line, U, and along the p coordinate line, V.” 


Contrarily in Ref. 4, it stated (for 3-D case) that 

“The contravariant velocity components U, V, W are in directions normal to constant 
£, r], C surfaces, respectively.” 

The above two statements mainly contradict each other the directions of contravariant 
velocity components. As we shall see, neither of these statements is strictly correct. The 
contravariant velocity components are neither the decomposed physical components of 
velocity vector nor in directions normal to constant £,r7,C surfaces, respectively. The 
objective of this article is to clarify what are the contravariant velocity components and 
to further explain their physical implications. 

To avoid confusion, we will start with defining terms based on the algebra of vectors. 
The term component of a vector will be defined. With the definition of covaraiant and 
contravariant base vector systems for a given coordinates, the general mathematical term 
“contravariant component will be discussed and its counter term covariant component 
will also be explained. For simplicity and clarity, we will use several 2-D linear cases as 
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examples. Here in the 2-D examples, the (a; 1 ,* 2 ) coordinates are interchangeable with 
(£,77) coordinates. 


II. Algebra of Vectors 

In an n-dimensional vector space a set of n linearly independent vectors e i? e 2) 

— e„ is called a base vector system. Any vector A in the space can be expressed as a unique 
combination of the base vectors A = c x ex + C2e2 H b c n e n , or using the summation 

convention, A = c*e*. The vectors c 1 e u c 2 e 2 , , c„e„ are component vectors of vector 

A in the e x ,e 2 , , and e„ directions, respectively. The coefficients, c x ,c 2 , » aild 

c„, are called components of vector A in the e x , e 2 , , and e n directions, respectively, 

with respect to the base (e x , e 2 , , e n ). 

For example, in Cartesian coordinates, a vector A = 5i + 7j implies 

e x = i, e 2 = j, ci = 5, c 2 = 7 

rewrite A = 1 0 (i — j)/2 + 12j, then that implies 

ei = (i - j)/2, e 2 = j, c x = 10, c 2 = 12 (2) 

rewrite A = — (2i) + 7 (i + j). then that implies 

ei = 2i, e 2 = i + j, c x = -1, c 2 = 7 (3) 

As illustrated above, the invariant vector A can be expressed in different base vector 
systems which may or may not be orthogonal and normalized. While the “component 
vectors” are physical vectors, the magnitude of the components of a vector not only depend 
on the directions but also depend on the magnitude of the base vector system. Obviously, 
without the prescription of “base vectors” , the “components’ of a vector have no meaning. 


III. Contravariant and Covariant Base Vectors 

Let x i be the coordinates of a point and r be a position vector. As shown in Fig. 1, 
a covariant base vector, g* = is tangent to its corresponding coordinate line, x\ 

and a contravariant base vector, g l = Vaff, is normal to the other coordinate lines x 3 
(or the surface of x i = constant), j ± i- For any given coordinates x l , its contravariant 
and covariant base vectors, g l and g, : are uniquely determined. They have the reciprocal 
relation, 


g l ■ Sj = 5 


h: 


if i ^ 3 
if i = 3 


(4) 


where is the Kronecker delta. Based on the reciprocal relation, Eq. (4), we have 
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g = 


gi I cos(gLgi) 


, (no summation in i). 


If they are orthogonal, then | g* |= j^-j and g l and gj are in the same direction. Note 
that, the sketches in Fig. 1 (and Fig. 2) only give a general indication of the directions of 
g i and g i and do not imply any thing about their magnitudes. 
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Actually, we have dv = g X dx^ and the differential of arc length ds is determined from 

ds 2 = dr • dr; therefore . 

ds 2 = g, ■ g j dx t dx j = g lJ dx l dx j 

Here, the quantities gij are components of a covariant tensor of rank two called the metric 
tensor. Similarly, We can define a contravariant tensor g i = g ■ gL Obviously, g is 
conjugate or reciprocal tensor of g^, i.e. 

9 lk 9k 3 = S) 

(g i • g j = 9 th gJb • gj = 9 tk 9kj = 5}) 

g lk and g^j are two fundamental tensors. Let g = \gij\ and G = \g l:l j are determinants of 
and. 5 then \J~Q -—l / \/G where \[Q is the Jctcobiciii of coordinate transformation from 
old to new and yfG is the Jacobian of coordinate transformation from new to old. 

IV. Contravariant and Covariant Components of a Vector 

A given vector A can be expressed in contravariant and covariant base vector systems 
as 

A = A'gi = Aig 1 (5) 

or 

A = A x gi + A 2 g 2 + A 3 g3 +A n gn (6) 

A = Aig 1 + A 2 g 2 + A 3 g 3 +A n g n (7) 

Here a super-script index is used for “contravariant” and a subscript index is used for 
“covariant”. According to the definitions in Sect. 1, A 1 gi,A 2 g 2 , ,A n g n are com- 
ponent vectors of vector A in the directions of gi,g 2 , >gm respectively. Simi- 
larly, Aig 1 , A 2 g 2 , , A n g n are component vectors of vector A in the directions of 

g^g 2 ^ ,g n 5 respectively. Accordingly, A 1 are components of vector A in the direc- 
tions of gj, with respect to the base vector system (gi,g 2 5 >gn)- Similarly, Ai are 

components of vector A in the directions of g 1 , with respect to the base vector system 

(g\g 2 > >g n )- As described in Sect. I, in general, the base vector system can be 

arbitrary (as long as the base system is complete). However, in Eq. (5), g » and g* are 
two uniquely defined base vector systems for a given coordinates x . gt is the covariant 
base vector system and g l is the contravariant base vector system for a given coordinates. 
Hence, we call 

A 1 are contravariant components of vector A in the directions of g t , with respect to 
the covariant base vector system (gi, g 2 , , g n)- 

A{ are covariant components of vector A in the directions of g*, with respect to the 
contravariant base vector system (g x ,g 2 , ,g n )- 

For simplicity, we call A 1 contravariant components of vector A and A x covariant 
components of vector A. 

When a vector is expressed in covariant and contravariant base vector systems, as shown in 
Eq. (5), its corresponding “components” are formed as contravariant and covariant vectors, 
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respectively. They follow the coordinate tr3.nsforin3.tion rule. Indeed, in Eqs. (6) and (7), 

A* = (A 1 , A 2 , , A n ) is a contravariant vector or a contravariant tensor of rank one, 

and A) = (Ai,A 2 ,~ ,A„ ) is a covariant vector, or a covariant tensor of rank one. 

Un- ambiguously, A^,A 2 - - - etc. are the components of the contravariant tensor A and 

Ai , A 2 etc. are the components of the covariant tensor Ai. Therefore, what we call 

“contravariant components” of a given vector are really components of a contravariant 
vector and, vis versa, what we call “covariant components” of a given vector are really 
components of a covariant vector. 

One can show that, from Eqs. (4) and (5), 

A 1 — A • g l and Ai = A • gi (8) 

Contravariant and covariant tensors, A*, A; can be transformed each other through the 
two fundamental tensors such as 

A 1 = g lj Aj, and Ai = gij A j 

Eq. (8) is useful for evaluation and commonly leads to a mis-interpretation that A is the 
component of vector A in the direction of g* and Ai the components of vector A in the 
direction of g*. This mis-interpretation stems from the common use of a dot-product with 
a unit vector, or “projection” , to calculate the component of a vector in the direction of 
that unit vector. (In projection, it implies that a vector is decomposed into two component 
vectors; one parallel to and another perpendicular to a given direction.) In Eq. (8), neither 
the contravariant base vectors nor the covariant base vectors are normalized, so that Eq. (8) 
is not a process of projection. Indeed, following the definition of Eq. (5). there should be 
no confusion and mis-interpretation as mentioned above. 


V. 2-D Examples and Continuity Equation 
Now in 2-D, let a position vector r(£, rj) with coordinate system (£, g), and 

A =- ui + v) = A 1 gi + A 2 g 2 = Aig 1 + A2g 2 
g 1 = = £ x i + j 

g 2 = Vg = g x \ + riyj 

dr . 

gl = ai = 1 + Vi J 

dr . 

g2 = ^ = Xr,l + y v J 

A 1 = (A • VO = £ x u + £ y v 
A 2 — (A • Vg) = g x u + g y v 
Ai = (A • gi) = x^u + y^v 
A 2 = (A • g 2 ) = x v u + y v v 
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Here A 1 is the contravariant component of A in the direction of £ and A 2 is the con- 
travariant component of A in the direction of 77, with respect to the base vector system 
\ a ig the covar i an t component of A in the direction of V£ and A 2 is the covariant 

component of A in the direction of V77, with respect to the base vector system (V£, V77). 

As one can see, neither contravariant nor covariant components are physical quanti- 
ties. Their relations to physical quantities can be illustrated, in 2-D, as Fig. 2. In Fig. 2, 
the vector A = OC is expressed by physical component vectors OB_ and OD _ in the direc- 
tions tangent to corresponding coordinate lines, gi, g 2 , or component vectors OE_ and OF_ 
in directions normal to the other coordinate lines, g 1 , g 2 , respectively. 

Note, in Fig. 2, OBCD is a parallelogram formed by BC_ || g2 and CD_ j| gi, 
and OECF is a parallelogram formed by EC || g 2 and CF || g 1 . By the definitions of 
contravariant and covariant base vectors, we have 

BC -L g\ CD. 1 g 2 , EC±g u CF±g 2 . 


Here, OG and OH are projections of vector OC in directions normal to the other coordinate 
lines, g 1 and g 2 , respectively. Note that, except in orthogonal coordinates, OG and OH 
are not “components” of vector OC . ( OG_ + OH 7^ OC , in general.) It is OG and GC or 
QE_ and HC_ that are component vectors of vector OC; i.e. OG + GC = OH_ + HQ = QC 
as shown in the sketch. 

Here we would like to illustrate the above discussions with several examples. For the 
example A, let a vector A = 5 i + 7 j, with coordinates (£,77) such as 

f = 2ar, V = V 


= £/2i + T7j 

g 1 = = £ x i + tyj = 2i, g 2 = Vr; = TJ x i + Tiyj=i 

dr . dr 

Si = af = 1/21 ’ S2 = s; =J 

Therefore, the contravariant components are 

A 1 = A • g 1 = 10, A 2 = A ■ g 2 = 7 
and the covariant components are 

A\ = A • gi = 5/2, A 2 = A • g 2 = 7. 


Fig.3 shows the sketch of this example, with magnitude and direction. Note that the new 
coordinates (£, 77) are same as the ( x , y) Cartesian coordinates, except £ is scaled by 2 x. 
While g 2 and g 2 are identical with in both magnitude and direction, gi and g 1 are in the 
same direction as x-coordinate, i, but different in magnitude, and the magnitudes of A 
and A\ are scaled accordingly. 

Example B, let the the same vector A = 5i + 7j, with new coordinate lines (Cv) so 
as 


£ = 2 x, rj = x + y 
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r(£,r?)=4/2i+(r?-£/2)j 
g 1 = 2 i, g 2 = (i+j) 
gi = (i — j)/2, g 2 =j 

Therefore, the contravariant components are 


A 1 = 10, A 2 = 12 

and the covariant components are 

A x = —1, A2 = 7. 

These are identical as Eqs. (2) and (3) in Sect. 1, written as 

A = A : gx + A 2 g 2 = 10(i - j)/2 + 12j 
A = Axg 1 + A 2 g 2 = — (2i) + 7(i + j) 

Fig. 4a shows detailed sketches of the base vectors, with magnitude and direction, and re- 
lations with ( x , y ). Fig. 4b shows the vector A expressed in its contravariant and covariant 
base vector systems. Since the magnitudes of base vectors are not normalized, the marked 
number are scaled to each base vectors accordingly. 

It is interesting to note that if vector A is decomposed into two component vectors 
so that one is parallel to g 1 and another is normal to g 1 , then A = 5/\/2 +7/2 

These are the results of projection on to unit vectors (i — j)/\/2and (i + j)/\/2 , and are 
quite different from the expression in contravariant or covariant base vector system. 

Example C, let a new coordinate system be £ = (x — y)/ 2, 77 = y, and hence 
r (£,77) = (2£ + 77)1 + 77J. Therefore, the contravariant and covariant base vectors are 

g 1 = (i -J)/ 2 , g 2 =j 

gx = 2 i, g 2 = (i+j) 

Hence, again as Eqs. (2) and (3), 

A = Axg 1 + A 2 g 2 = 10(i — j)/2 + 12 j 


A = A 1 g 1 + A 2 g 2 = — (2i) + 7 (i + j) 


Fig. 5 shows the contravariant and covariant base vectors for this new coordinate system. 
Note that the coordinate system in Example C is the conjugate or reciprocal of the the co- 
ordinate system in Example B. The contravariant base vectors in Example C are covariant 
base vectors in Example 13 and the covariant base vectors in Example C are contravariant 
base vectors in Example B. Clearly, the vector A is an invariant and the specification of 
coordinate system is essential for vector decomposition. 


In fluid dynamics, the continuity equation can be written in integral and differential 
forms as 

/ p dv + / pW-ndS = 0, 

dt J v Js 
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(9) 


fiPiff)* + ^PW • (g)l g 1 + ^pW • (g)i g 2 = 0, 

^p(p)^ + ^p(9)*u + -T^pig^v = o. 

Here, p is density, W = ui + vj the velocity vector, and (7 = W • g 1 and V = W • g 2 are 
contravariant velocity components in the directions of covariant base vectors gi and g2, 
respectively. Clearly, here ^/g is the same as J -1 in Eq. (1), the Jacobian of transformation 

from old to new coordinates (or yfg~ l — VG = J). In a finite volume formation, g is the 
volume of the grid cell and (g)? g 1 and (g)? g 2 are surface normal in the directions of g 1 
and g 2 , respectively. The concept of mass conservation can be interpreted accordingly. As 
Eq. (9), while the contravariant velocity ( U , V) is closely associated with the dot-product of 
velocity with surface normal, there should not be confusion that the contravariant velocity 
components, U, V, are velocity components in the directions of normal to the corresponding 
surfaces. 

Finally, it is important to take note that the contravariant velocity components, 
(17, V), of the the velocity vector W do not have the physical dimension of W (velocity) 
itself. In a discretized numerical grid domain, contravariant velocity components, (U,V), 
have dimension of velocity divided by a length scale. Conversely, covariant velocity compo- 
nents have dimension of velocity multiplied by a length scale. They magnitudes are scaled 
by their corresponding base vectors accordingly. 

In fact, in addition to their magnitudes, the dimensions of contravariant and covari- 
ant components are controlled by the dimensions of the covariant and contravariant base 
vectors, respectively. In Eq. (4) the contravariant and covariant base vectors are reciprocal 
not only in magnitude but also in dimension. This can be clearly illustrated by writing 
continuity equation in the cylindrical coordinates as 


£ = r = y/(a: 2 + y 2 ), rj = 0 = tan l x/y 


gi ~ sinO i + cosd j, g 2 = r cosd i — r sin9 j 
g 1 = sinO i + cosd j , g 2 = 1/r cosB i — \/r sin6 j 
U = u sinO + v cosP, V — 1 /r(u cosQ — v sind) 
cf c) d 

— p/r + —p/r(u sinO + v cosO ) -I- wt:P/c 2 (u cosO — v sinO) = 0. 
ot or oO 

Obviously, in the above, the dimension of contravariant velocity component U is different 
from that of contravariant velocity component V. 

VI. Concluding Remarks 


In this paper we have reviewed the basics of tensor analysis in an attempt to clarify 
some misconception regard contravariant and covariant vector components as used in fluid 
dynamics. We have indicated that contravariant components are components of a given 
vector expressed as a unique combination of the covariant base vector system and vis versa, 
the covariant components are components of a vector expressed with the contravariant base 
vector system. Mathematically, expressing a vector with a combination of base vector is a 
decomposition process for a specific base vector system. Hence, the contravariant velocity 
components are decomposed components of velocity vector along the directions of coordi- 
nate lines, with respect to the covariant base vector system. However, the contravariant 
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(and covariant) components are not physical quantities. Their magnitudes and dimensions 
are controlled by their corresponding covariant (and contravariant) base vectors. 
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direction only). 
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J 


A = PC = OB + OP = OE + OF 
OB = (A • g J )gi = A 1 gi 
OD = (A • g 2 )g 2 - A 2 g 2 
OE = (A • g^g 1 = Ai g 1 
QF = (A • g 2 )g 2 = A 2 g 2 

0G = ( A.gWhWl 

OF=(A-g 2 )/|g 2 | = A 2 /|g 2 | 

Fig. 2 The relations among its physical quantities of a given vector A to contravariant and 
covariant components 
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%- Z 4- 6 § )° 

£ = 2x, r) = y 

?) = f/2i 4- ry 

g 1 = v? = y + y = 2i, g 2 = Vt, = 17,1 + y = j 
«> = = V2i, g2 = |^ = j, v 5 = l/ 2 

Fig. 3 The Contravariant and Covariant Base Vectors in the Coordinates of £ = 2x and 
T] = y, with Vector A = 5i + 7j 
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? = 2x, 7 ? = x + y 

r(£,*7) = £/2i+fo-e/2)j 
g 1 — 2 i, g 2 = (i+j) 

Si = (i ~ j)/2, g2 = j, = 1/2 


\ 

J 



Fig. 4a The Contravariant and Covariant Base Vectors in the Coordinates of £ = 2x and 




* 


A == 10 gl + 12g 2 = 10(1 - j)/2 + 12j J* 

A = — g 1 4- 7g 2 = — (2i) + 7(i + j) 

Fig. 4b A Vector A = 5i + 7j Expressed with respect to Contravariant and Covariant Base 
Vector Systems in the Coordinates of £ = 2x and r? = x + y 
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